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A weak-coupling formalism for superconducting states possessing both singlet (even parity) and 
triplet (odd parity) components of the order parameter in centrosymmetric crystals is developed. It 
is shown that the quasiparticle energy spectrum may be non-degenerate even if the triplet component 
is unitary. The superconducting gap of a mixed-parity state may have line nodes in the strong spin- 
orbit coupling limit. The pseudospin carried by the superconducting electrons is calculated, from 
which follows a prediction of a kink anomaly in the temperature dependence of muon spin relaxation 
rate. The anomaly occurs at the phase boundary between the bare triplet and mixed-parity states. 
The stability of mixed-parity states is discussed within Ginzburg-Landau theory. The results may 
have immediate application to the superconducting series Pr(Osi_a;Rua;)4Sbi2. 
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I. INTRODUCTION 

Very recently, magnetic susceptibility and electrical 
resistivity measurements on the filled skutterudite se- 
ries Pr(Os 1 _ a .Ru a; )4Sbi2 found that superconductivity 
exists in the whole concentration range < x < li 
While PrRu4Sbi2 is a conventional s-wave superconduc- 
tor with T c =l K, PrOs4Sbi2 is a heavy fermion ma- 
terial 3 with X" c =1.85 K. There is experimental evidence 
that PrOs4Sbi2 has two superconducting (SC) phase d 3 ! 4 ' 5 
and that the quasiparticle energy spectrum has point 
nodes m& These facts indicate that the SC order param- 
eter is definitely unconventional. Zero-field muon-spin 
relaxation (/zSR) measurements suggest that the pairing 
in PrOs4Sbi2 can be triplet (odd parity)^ The compe- 
tition between the conventional s-wave and triplet order 
parameters may result in the appearance of a new SC 
state in Pi^Osi-^Ru^Sb^ which has both singlet and 
triplet components and therefore is a mixed-parity (MP) 
state. 

So far, several systems in which MP SC states may 
occur have been studied theoretically. Mineev and 
Samokhin^ showed that the normal state may be unsta- 
ble with respect to a helical MP structure if the product 
of the two different parity representations of the symme- 
try group contains a vector representation. The possi- 
ble occurrence of MP states due to inhomogeneity of the 
order parameter in the Larkin-Ovchinnikov-Fulde-Ferrel 
state was addressed in Refs. IflUlfl MP states in two- 
dimensional superconductors, in which inversion symme- 
try is broken due to low dimensionality, were studied by 
Gor'kov and Rashba. 11 Finally, an intriguing prediction 
was made by Volkov et alM^ about the possibility of 
generating a triplet condensate in mesoscopic ferromag- 
net/singlet superconductor multilayers. 

In this paper, we analyse MP SC states in bulk crys- 
tals with inversion symmetry in zero magnetic field and 
describe the possible experimental manifestation of the 
MP state in P^Osi-^Ru^Sb^. We discuss the pos- 
sibility of the realization of a MP state in terms of 



Ginzburg-Landau (GL) theory. We develop a weak- 
coupling formlaism for MP states based on the general- 
ized BCS approach We obtain the quasiparticle energy 
spectrum and show that it may be non-degenerate even 
if the triplet component is unitary. It is shown that the 
gap function of the MP state in the strong spin-orbit 
limit may have line nodes, in contrast to bare triplet 
states. The Gor'kov formalism is used to obtain the 
self-consistent equation of the gap function and to calcu- 
late the pseudospin carried by the SC electrons. A kink 
anomaly in the temperature dependence of the muon spin 
relaxation rate is expected on the boundary between the 
bare triplet and MP states. 

II. GINSBURG-LANDAU MODEL 

The explicit form of the GL functional can be estab- 
lished only if the transformation properties of the multi- 
dimensional order parameters with respect to the spa- 
tial symmetry transformations are known. Thus far, 
the symmetry of the order parameter in PrOs4Sbi2 has 
not been determined unambiguously. 15 Also, the effects 
of multi-dimensionality of the order parameter are be- 
yond the scope of this paper. Thus, we adopt an effec- 
tive one-component model which only takes into account 
gauge, inversion and time reversal but not the crystallo- 
graphic symmetry. We denote by n = l^e 1 ^ 1 the non- 
vanishing component of the singlet order parameter and 
by £ = |£|e J ^ 2 that of the triplet order parameter. A ho- 
mogeneous situation in zero magnetic field is considered 
in the following, so that the gradient terms of n and £ 
are neglected. The GL potential is 

F = ai [T -T cl (x)}\ V \ 2 +h[T -T c2 (x)]\e (I) 

+a 2 \v\ 4 + M£| 4 + ciM 2 |£| 2 + c 2 ( V 2 C 2 + V* 2 e), 

where T is temperature and T c \ t2 (x) are concentration 
dependent transition temperatures for the singlet and 
triplet order parameters, respectively. The latter can be 
roughly fitted to the experimental datai as linear func- 
tions T cl (x) =T s [l-a(l-x)], T c2 (x) =T t (l-j3x) with 
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T s = 1.2 K, T t = 1.8 K, a = 1.07, and (3 = 1.03. The rest 
of the parameters in are assumed to be undetermined 
constants^ 

F can be immediately minimized with respect to the 
relative phase A<p = 4>i — 4> 2 . Depending on the sign of 
c 2 , the following MP states are energetically favoured 

= 71-/2,371-/2, r)£* =±i|?7||e| for c 2 > 0; (2a) 
A0 = O,tt, rit* = ±\ri\\£\ forc 2 <0. (2b) 

The phase diagram of adapted for 

Pr(Osi_ K Ru a; )4Sbi2 is shown in Fig. If 
4a 2 6 2 > (ci - 2|c 2 |) 2 , one of the MP states © ex- 
ists in the sector bounded by two second-order transition 
lines 

caiT cl (x) - 2a 2 (5iT c2 {x) 



T{x) = 
on the singlet side and 



cai — 2et 2 /3i 



T(x) 



2b 2 a 1 T cl (x) - c/3iT c2 {x) 
2b 2 a± — c[3\ 



(3) 



(4) 



on the triplet side. If 4a 2 b 2 < (c x - 2|c 2 |) 2 , the MP 
states (01 do not minimize F and a first-order phase tran- 
sition between the bare singlet and triplet states occurs 
at the line 



T(x) = 



aiVb~2~T c x(x) - f3 ly /a^T c2 (x) 



(5) 



Phenomenologically, therefore, both types of the phase 
diagram shown in Fig. ^ may occur with the same 
probability in the sense that they are described by the 
same fourth-order model QJ. Another doped system 
Ui-^Th^Beis is an example in which the co-existcncc 
of two SC order parameters is preferred to the first-order 
phase transition scenario JLi£ 



III. GENERALIZED BCS APPROACH 

A realistic description of superconductivity in 
Pr(Osi_ 2; Ru a ;)4Sbi2 should be based on a strong- 
coupling approach which would incorporate the super- 
conductivity mechanism and coupling of SC electrons 
to the Os(Ru) ions. However, at present even the su- 
perconductivity mechanism of pure PrOs4Sbi 2 remains 
vague. In this paper our goal is to elucidate the quali- 
tative features of MP states rather than to make quan- 
titative predictions. Hence, we use the weak-coupling 
approach, implicitly assuming that the main ingredients 
of the theory, the band energy relative to the chemical 
potential e(k) and the pairing interaction matrix element 
V S1 . 



S4 (k, k'), depend on the Ru concentration x. 
The generalized mean-field BCS Hamiltonian i: 



,14 
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k,Sl,S2 



- A s 1S2 (- k )a-k Sl a kS2 ], 



FIG. 1: Sketch of the phase diagram of Pr(Osi_ a; Ru :E )4Sbi2. 
Solid lines represent a fit to the experimental datai The other 
lines are calculated from the model Q. (a) Broken lines are 
the boundaries of the MP state for 4a2&2 > (ci — 2 1 C2 1 ) 2 [See 
Eqs. 10 and Q)]. (b) Dash-dot line represents a first-order 
phase boundary for 4a2&2 < (ci — 2|c2 1) 2 [Eq. ©]. 
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and 



where aj^ (ak s ) is the creation (annihilation) operator of 
an electron with wave vector k and pseudospi] 
the gap function matrix is 

A SlS2 (k) = 



V S2SlS3Si (k,k')(a k , S3 a- k , Si ). (7) 

k',s 3 ,s 4 



The gap function is conveniently parametrized by an even 
scalar function ip(k) = r]iipi(k) for the singlet compo- 
nent and an odd vectorial function d(k) = ^j&cL^k) 
for the triplet component. Here ipi(k) and dj(k) are 
basis functions of the irreducible representations of the 
point group and rji and & are the corresponding order 
parameters^ In the MP state the gap function takes 
the form 



A(k) = [^(kjffo + d(k)a]ia y , 



(8) 



where <7o is the 2x2 unit matrix and <x = (a x ,a y ,az) 
are the Pauli matrices. 

The Bogoliubov transformation for the diagonalization 
of © yields the quasiparticle energy spectrum 



(9) 



E k± = v/e(k) 2 + A±(k) 2 , 

with two gaps A± (k) defined by 

A ± (k) 2 = |^(k)| 2 + |d(k)| 2 ± |p(k) + q(k)|, (10) 

where two real vectors p(k) = ^(k)d*(k) + -0*(k)d(k) 
and q(k) = i[d(k) x d*(k)] are introduced. They are 
orthogonal to each other, therefore 



|p(k)+q(k)| = v/p(k) 2 +q(k) 



(11) 



(6) 



Note that if q(k) ^ 0, i. e. d(k) and d*(k) are not 
collinear, then p(k) is finite in the MP state. 

It follows that the quasiparticle spectrum JjJJ is degen- 
erate only if p(k) = q(k) = 0. Such a state corresponds 
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FIG. 2: (Colour online) Line nodes in the SC gap of the time- 
reversal invariant admixture of s- and p-wave states l|lMfl for 

161 > 161- (a) \v\ < 161; (b) \r,\ = (c) 161 < M < 161- 




to l|2a|) . It is not invariant with respect to time rever- 
sal K,, but it possesses a combined symmetry element 
U\{-k)IK,, where U\ is the gauge transformation and / is 
inversion. In a time reversal invariant MP state, q(k) = 
andp(k) = 2|V>(k)||d(k)| [See J20]. Then, the gap ac- 
quires the form 

A ± (k) 2 = (|V(k)| ± |d(k)|) 2 . (12) 

As follows from l|10|l . the gap of a MP state vanishes 
if both singlet and triplet components have nodes, i. e. 
tp(k) = and |d(k)| 2 = |d(k) x d*(k)|. In the strong 
spin-orbit coupling limitj2& these nodes can be found 
only at isolated points on the Fermi surfac o 22 i 23 unless 
the Fermi surface crosses the boundary of the Brillouin 
zone^ However, A_ may have additional line nodes if 
p 7^ 0, as illustrated by the following example. 

Keeping in mind the possible application to 
Pr(Osi_ K Ru x )4Sbi2, we consider a p-wave state 
which may arise under Th crystallographic symmetry. 
Time reversal symmetry is broken in the SC state of 
PrOs4Sbi2i£ Here, for the sake of simplicity, a time- 
reversal invariant admixture of s- and p-wave states is 
considered, in which the gap takes the simple form (11211 . 
We can show that this simplification does not affect the 
main result, line nodes may exist for finite q as well. 

Let us take a p-wave state belonging to the irreducible 
representation T u of Th- This representation has two 
linearly independent sets of basis functions^ Therefore, 
in general, one should take into account two order pa- 
rameters of the same symmetry. We consider the state 
denoted (0, 0, 1) in the order parameter spaced The MP 



gap function is given by 

V>(k) = 7 h d(k) = £ifc x y + iih y % (13) 

where n, £1 and £2 are the T and x dependent order 
parameters. A_ (k) vanishes at two lines defined by the 
equation 

\ti\ 2 kl + \&\ 2 kl = \il\ 2 - (14) 

The line nodes exist while |7/| < max(|£i |, I^D- This 
condition is obeyed close to the boundary with the bare 
triplet state. Fig.|2shows the line nodes on a unit spher- 
ical Fermi surface for |6| > 1^2 1- Note that in this ex- 
ample the line nodes are not remnants of line nodes of 
the singlet state, which is fully gapped, but rather they 
appear due to a distortion of point nodes of the triplet 
component. 



IV. GREEN'S FUNCTIONS AND GAP 
EQUATION 

Beginning with the equations of motion for Heisen- 
berg operators, the Gor'kov equations 2 ^ for the normal 
G ss /(k, r) and anomalous F ss /(k, r), Fj s , (k, r) tempera- 
ture Green's functions in the MP state are obtained in 
the usual wayi2S^£ The Green's functions are defined as 

G^(k,r) = -(T T {a ks (r)al,(0)}) 
F ss ,{k,r) = (T r {a ks (r) a _ ks ,(0)}) 
Fj s ,(k,r) = <T T { a t_ ks(r )4 s ,(0)}), (15) 

where T T is the imaginary time ordering operator. The 
Fourier transform is A(k, r) = A(k, w„)e _lu,, ' r , 

where A stands for G, F, or F'; T is temperature and 
u> n = ttT(2ii+1) is the Matsubara frequency for fermions. 

The resulting Gor'kov equations for the Fourier trans- 
forms are 

[^„-e(k)]G(k, Wn ) + A(k)F t (k,c i ; n ) = a , 
[iw„ + £ (k)]Ft(k,a; n ) + At(k)G(k,a; n ) = 0. (16) 

They are solved by 



G(k,w„) 
F(k,w„) 



e(k) 2 + |^(k)| 2 + |d(k)| 2 ]i?o + [p(k) + q(k)]<7 . 



:(k) 2 



l^(k)| 2 



d(k)| 2 ]A(k) - 



<k)], 



|p(k)+q(k)|fi(k) 



) 



(17) 



where The denominator introduced in the definition of f2(k) 

in 1 n \^ 1 r / /i \ n \ 1 ■ n \ „ j/i ni allows one to put the self-consistent gap equation in a 

n, n d(k)p(k)(To + V(k)PW +«q(k) x d(k) <r,_ . v BV H 

m = | P (k) + q(k)| M 

(18) 
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concise form. Using (J7J), l|15(l . and (|17() one obtains 
A slS2 (k) = - J2 V S2S1S3S4 (k,k')T S3S4 (k',T), (19) 

k',s 3 ,s 4 



where the indices and JJ- denote the pseudospin-up and 
down states, respectively. Using Eqs. fT5|l . (fTTjl . J22J), and 

the oddness of p(k), one obtains the average pseudospin 
of the Cooper pairs, 



where 



HKT) = ^±«tanh^ k 



4£ k 



2T 



A(k) - fi(k) , /£ k _ . 

; tanh ( ) . (20) 



4£ k 



2T 



s(t)=x; 



e(k)q(k) 
p(k)+q(k)| 



tanh(£; k _/2T) 

tanh(£ k+ /2T) 
4Ek7 



(23) 



V. MUON SPIN ROTATION IN 
MIXED-PARITY STATES 

fiSK measurements are widely used to reveal the na- 
ture of a SC state as they provide invaluable informa- 
tion about the distribution of the local magnetic field H. 
There are several possible sources of internal magnetic 
fields in superconductors. They include magnetic mo- 
ments of localized states, two-dimensional imperfections 
such as domain walls and the surfaceji^2£ and magnetic 
moment of Cooper pairs. Sigrist and RiceA^ showed that 
the latter contribution always vanishes in singlet pairing 
states, and it is non-vanishing in bare triplet states only 
if q(k) 7^ 0. The zero- field /iSR spectra are usually fitted 
using of the Kubo-Toyabe function^ 

9(t) = 1 + 1(1 -5 2 t 2 )exp(-i<5 2 t 2 ), (21) 

where t is time and the relaxation rate 5 is proportional 
to J (H 2 ). The temperature dependence of 6 in the MP 
states can be obtained using the Green's functions (|T7jl . 

A calculation of the actual magnetic moment with 
strong spin-orbit coupling requires detailed knowledge 
of the electronic band structure. At every point k in 
momentum space, the operators a k and a ks can be ex- 
pressed as linear combinations of Bloch sums consisting 
of orbital and spin parts. In the calculation of the mag- 
netic moment, a quadratic form of the coefficients of this 
transformation will be an additional factor in the sum- 
mand of the resulting expression l)23|l (see below). How- 
ever, the qualitative features of S(T) may be discerned by 
calculating the average pseudospin of the Cooper pairs^S 
The pseudospin operator is defined as 

k 
k 
k 



Hence, as in the case of a bare triplet state, S(T) van- 
ishes if q(k) = 0. If the bare triplet state is nonunitary 
[q(k) 7^ 0] then, as follows from (|23|l , the temperature de- 
pendence of 5 has a kink at the point of the bare triplet- 
to-MP phase transition with a discontinuity in dS/dT 
proportional to <ip(k) 2 /eff. 



VI. SUMMARY 



Using a GL-type model, we showed that a MP SC state 
can be stable in a finite region of the T-x phase diagram 
of Pr(Osi_ x Ru :c )4Sbi2. The phase transitions from the 
bare singlet and triplet states to the MP state are second- 
order. The general formalism describing the MP states 
within BCS approach is developed. The quasiparticle 
energy spectrum is degenerate only if p(k) = q(k) = 0. 
The gap in the MP state can have line nodes in strong 
spin-orbit coupling limit if p(k) 0. Normal and anoma- 
lous Green's functions and the self-consistent gap equa- 
tion are obtained. The MP state may be experimentally 
indicated by line nodes in the gap of the quasiparticle 
spectrum and a kink in /xSR relaxation rate S(T). 
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